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Abstract. We consider i.i.d. random variables {oj(b) : b G E^} parameterized by the 
family of bonds in Z d , d ^ 2. The random variable w(b) is thought of as the conductance 
of bond b and it ranges in a finite interval [0, Co] . Assuming the probability m of the event 
{uj(b) > 0} to be supercritical and denoting by C(w) the unique infinite cluster associated 
to the bonds with positive conductance, we study the zero range process on C(co) with 
ct;(6)-proportional probability rate of jumps along bond b. For almost all realizations of 
the environment we prove that the hydrodynamic behavior of the zero range process is 
governed by the nonlinear heat equation dtp = mV ■ (W(/>(p/m)), where the matrix T> 
and the function (j> are oj-independent. As byproduct of the above result and the blocking 
effect of the finite clusters, we discuss the bulk behavior of the zero range process on Z d 
with conductance field uj. We do not require any ellipticity condition. 

Key words: disordered system, bond percolation, zero range process, hydrodynamic limit, 

homogenization, stochastic domination. 
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1. Introduction 

Percolation provides a simple and, at the same time, very rich model of disordered 
medium [G], |Kej . The motion of a random walker on percolation clusters has been 
deeply investigated in Physics (see |BH] and reference therein) and also numerous rigorous 
results are now available. In the last years, for the supercritical percolation cluster it has 
been possible to prove the convergence of the diffusively rescaled random walk to the 
Brownian motion for almost all realizations of the percolation |SS| , [BBj , [MPJ , improving 
the annealed invariance principle obtained in [DFGWj . We address here our attention to 
interacting random walkers, moving on the supercritical Bernoulli bond percolation cluster 
with additional environmental disorder given by random conductances (for recent results 
on random walks among random conductances see |BPj . [M] , [F] and reference therein). 

Particle interactions can be of different kind. An example is given by site exclusion, the 
hydrodynamic behavior of the resulting exclusion process has been studied in [F]. Another 
basic example, considered here, is the zero range interaction: particles lie on the sites of 
the infinite cluster without any constraint, while the probability rate of the jump of a 
particle from site x to a neighboring site y is given by g(rj(x))uj(x, y), where g is a suitable 
function on N, rj(x) is the number of particles at site x and co(x,y) is the conductance 
of the bond {x,y}. We suppose that the concuctances are i.i.d. random variables taking 
value in [0, cq}. 

The above exclusion and zero range processes are non-gradient systems, since due to 
the disorder the algebraic local current cannot be written as spatial gradient of some local 
function. Nevertheless, thanks to the independence of the conductances from any bond 
orientation, one can study the hydrodynamic behavior avoiding the heavy machinery of 
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non-gradient particle systems [V], [KL] [Chapter VII]. Indeed, in the case of exclusion pro- 
cesses, due to the above symmetry of the conductance field the infinitesimal variation of 
the occupancy number rj{x) is a linear combination of occupancy numbers. This degree 
conservation strongly simplifies the analysis of the limiting behavior of the random empir- 



processes, this degree conservation is broken. Nevertheless, due to the symmetry of the 
conductance field, adapting the method of the corrected empirical measure [GJlj to the 
present contest one can reduce the proof of the hydrodynamic limit to an homogenization 
problem plus the proof of the Replacement Lemma. The resulting diffusive hydrody- 
namic equation does not depend on the environment and keeps memory on the particle 
interaction. 

The homogenization problem has been solved in [F] also for more general random con- 
ductance fields. The core of the problem here is the proof of the Replacement Lemma. 
This technical lemma compares the particle density on microscopic boxes with the parti- 
cle density on macroscopic boxes and it is a key tool in order to go from the microscopic 
scale to the macroscopic one. This comparison is usually made by moving particles along 
macroscopic paths by microscopic steps and then summing the local variations at each 
step. The resulting method corresponds to the so called Moving Particle Lemma and be- 
comes efficient if the chosen macroscopic paths allow a spread-out particle flux, without 
any concentration in some special bond. While for a. a. oj any two points x, y in a box A^v 
of side N centered at the origin can be connected inside the infinite cluster by a path 7 X) „ 
of length at most O(N) [2E], it is very hard (maybe impossible) to exhibit such a family of 
paths { r ix,y}x,yeK N with a reasonable upper bound of the number of paths going through 
a given bond b, uniformly in b. Due to this obstacle, we will prove the Moving Particle 
Lemma not in its standard form, but in a weaker form, allowing anyway to complete the 
proof of the Replacement Lemma. We point out that in this step we use some technical 
results of |APj . where the chemical distance inside the supercritical Bernoulli bond perco- 
lation cluster is studied. It is only here that we need the hypothesis of i.i.d. conductances. 
Extending part of the results of [AP] , one would get the hydrodynamic limit of zero range 
processes among random conductances on infinite clusters of more general conductance 
fields as in [F], 

We comment another technical problem we had to handle with. The discussion in [GJlj 
refers to the zero range process on a finite toroidal grid with conductances bounded above 
and below by some positive constants, and some steps cannot work here due to the presence 
of infinite particles. A particular care has to be devoted to the control of phenomena of 
particle concentration and slightly stronger homogenization results are required. 

Finally, in the Appendix we discuss the bulk behavior of the zero range process on Z rf 
with i.i.d. random conductances in [0, cq], in the case of initial distributions with slowly 
varying parameter. Due to the blocking effect of the clusters with finite size, the bulk 
behavior is not described by a nonlinear heat equation. 

We recall that the problem of density fluctuations for the zero range process on the su- 
percritical Bernoulli bond percolation cluster with constant conductances has been studied 
in [GJ2]. Recently, the hydrodynamic limit of other interacting particle systems on Z d , 
or fractal spaces, with random conductances has been proved (cf. [FT], |F JL| . [JL| . |LF| . 
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[Val] ). We point out the pioneering paper [FrJ, where J. Fritz has proved the hydrody- 
namic behaviour of a one-dimensional Ginzburg-Landau model with conservation law in 
the presence of random conductances. 

2. Models and results 

2.1. The environment. The environment co modeling the disordered medium is given by 
i.i.d. random variables \oj(b) : b £ E^j , parameterized by the set E^ of non-oriented bonds 
in Z d , d ^ 2. to and co(b) are thought of as the conductance field and the conductance at 
bond b, respectively. We call Q the law of the field ui and we assume that co(b) G [0, cq] for 
Q-a.a. co, for some fixed positive constant cq. Hence, without loss of generality, we can 
suppose that Q is a probability measure on the product space $7 := [0,co] Ed . Moreover, 
in order to simplify the notation, we write co(x,y) for the conductance cj(b) if b = {x,y}. 
Note that co(x, y) = co(y, x). 

Consider the random graph G(co) = (V(co), E(co)) with vertex set V(co) and bond set 
E(co) defined as 

E(u) := {b£K d : co(b) > 0} , 

V(co) := {x £ Z d : x £ b for some b £ E(co)}. 

Assuming the probability Q(co(b) > 0) to be supercritical, the translation invariant Borel 
subset Qq C given by the configurations co for which the graph G(co) has a unique infinite 
connected component (cluster) C(co) C V(co) has Q-probability 1 [G] . Below, we denote 
by £{cj) the bonds in E{co) connecting points of C(co) and we will often understand the 
fact that co £ Qq. 

For later use, given c > we define the random field co c = (co c (b) : b £ E^) as 



otherwise . 



For c = we simply set co := coq. 



2.2. The zero range process on the infinite cluster C(co). We fix a nondecreasing 
function g : N — )■ [0, oo) such that ^(0) = 0, g(k) > for all k > and 

g* :=snp\g(k + l)-g(k)\<oo. (2.2) 

fceN 

Given a realization co of the environment, we consider the zero range process r\t on the 
graph G(co) = (C(co),£(co)) where a particle jumps from x to y with rate g(r](x))co(x,y). 
This is the Markov process with paths r](t) in the Skohorod space D([0,oo),N c ^) whose 
Markov generator C acts on local functions as 

£/fa) = E E 9(v(x))oo(x, x + e)(/(^ +e ) -/(r?)) , (2.3) 

e£B x£C(uj): 

where B = {±ei, ±e2, . . . , ie^}, ei, . . . , e<j being the canonical basis of Z d , and where in 
general 

rj(x) — 1, if z = x , 
jf'^jg) = r/(y) + 1, if z = y, 
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We recall that a function / is called local if f(rj) depends only on r](x) for a finite number 
of sites x. Since C(u) is infinite, the above process is well defined only for suitable initial 
distribution. As discussed in PQ, the process is well defined when the initial distribution 
has support on configurations rj such that \\rj\\ := ^2xeC(u) v( x ) a ( x ) < 00 > °(") being a 
strictly positive real valued function on C{oS) such that 

lo(x, x + e)a(x + e) ^ Ma(x) 

x£C{uj) e£B 

for some positive constant M. 

Given <p > 0, set Z(<p) :=Ek>o vVsW wh ere 5 (0)! = 1, ff(fc)! = ^(1)^(2) • • • 5 (fc) for 
fc ^ 1. Since Z(<^) is an increasing function and g(k)\ #(l) fe , there exists a critical value 
(p c £ (0, oo] such that Z(<p) < oo if (p < ip c and Z(jp) = oo if > </? c . Then, for ^ (p < cp c 
we define as the product probability measure on N c(tj) such that 

1 ip k 

v v (i){x) = k)= , fceN, xeC(u). 

Z{<p)g{k)\ 

Taking for example ol(x) = e - ^' in the definition of \\r]\\ one obtains that ^(H^H) < oo, 
thus implying that the zero range process is well defined whenever the initial distribution 
is given by u v or by a probability measure fj, stochastically dominated by v v . In this last 
case, as proven in [X], by the monotonicity of g one obtains that the zero range process 
r\t starting from \x is stochastically dominated by the zero range process Q starting from 
Up, i.e. one can construct on an enlarged probability space both processes rjt and Q 
s.t. rjt(x) ^ Ct(x) almost surely. Finally, we recall that all measures are reversible 
for the zero range process and that Ptp(e^ x ') < oo for some 6 > 0, thus implying that 
^(v(x) k ) < oo for all k > (cf. Section 2.3 of [KL] 1. 

As usually done, we assume that lim (/p ^ c Z(ip) = oo. Then, cf. Section 2.3 in [KL], the 
function R(np) := ^(77(0)) is strictly increasing and gives a bijection from [0, ip c ) to [0, 00). 
Given p S [0, 00) we will write (p(p) for the unique value such that R(<p) = p- Then we set 

v p '•= v<p(j>) , 4>(p) -■= vp(g{vx)) xeC(u). (2.4) 

As proven in Section 2.3 of |KL| . cj> is Lipschitz with constant g* . 

2.3. The hydrodynamic limit. Given an integer N ^ 1 and a probability measure p N 
on N^), we denote by 1P W the law of the zero range process with generator N 2 C (see 
(|2.3p ) and with initial distribution pL (assuming this dynamics to be admissible). We 
denote by E^ ^jv the associated expectation. In order to state the hydrodynamic limit, 
we define B(£l) as the family of bounded Borel functions on f2 and let T> be the d x d 
symmetric matrix characterized by the variational formula 

{a,Va) = - inf J V / w(0, e)(o e + V(r e w) - VM)\ eeCH Q(du;) \ , Va£R d , 

(2-5) 

where denotes the canonical basis of 

m:=Q(0£C(w)) (2.6) 

and the translated environment r e oj is defined as T e co(x, y) = u{x + e, y + e) for all bonds 
{x,y} in Erf. In general, denotes the characteristic function of A. 



ZERO RANGE PROCESSES ON THE SUPERCRITICAL PERCOLATION CLUSTER 



5 



The above matrix T> is the diffusion matrix of the random walk among random conduc- 
tances on the supercritical percolation cluster and it equals the identity matrix multiplied 
by a positive constant (see the discussion in [F] and references therein). 

Theorem 2.1. For Q-almost all environments u the following holds. Let po : R rf — > [0, oo) 
be a bounded Borel function and let {pn}n ^ i be a sequence of probability measures on 
j^C(lj) suc /j j or a ll 5 > and all continuous functions G on R rf with compact support 
(shortly G £C c (R d )), it holds 



lim fi N 

jVfoo 



N~ d V G(x/N)r](x)- [ G(x)p (x)dx 
77! * JR d 



><5 =0. (2.7) 



x&C(ui) 



Moreover, suppose that there exist po,p*,Co > such that p N is stochastically dominated 
by v po and the entropy H{p N \v P:t ) is bounded by CoN d . 
Then, for all t > 0, G £ C c (R d ) and 5 > 0, it holds 



lim P w N 



>6 =0, (2.8) 



N~ d V G(x/N)r] t (x) - [ G(x)p(x,t)dx 
771 \ 

where p : M. d x [0, oo) — >• R is assumed to be the unique weak solution of the heat equation 

d t p = mV • {VV4>(p/m)) (2.9) 
with boundary condition po at t = 0. 

We define the empirical measure n (rj) associated to the particle configuration rj as 



» I I I I I I I --. ■ i. — .■■ >/ I I I ■■■ 

N d 



where A^(R ) denotes the Polish space of non-negative Radon measures on R endowed 
with the vague topology (namely, v n — > v in Ai if and only if f n (f) — > v{f) for each 
/ G C c (R d )). We refer to the Appendix of [SJ for a detailed discussion about the space M. 
endowed of the vague topology. We write 71"^ for the empirical measure ir N (%), f]t being 
the zero range process with generator N 2 C. Then condition (|2.7p simply means that under 
p N the random measure ir N converges in probability to p^{x)dx, while under jv the 
random measure ir^ converges in probability to p(x, t)dx, for each fixed t ^ 0. In order to 
prove the conclusion of Theorem 12.11 one only needs to show that the law of the random 
path tt n £ D([Q,T},M) weakly converges to the delta distribution concentrated on the 
path [0, T] B t — > p(x,t)dx £ M (see |KL| [Chapter 5]). It is this stronger result that we 
prove here. 

Let us give some comments on our assumptions. We have restricted to increasing func- 
tions g in order to assure attractiveness and therefore that the dynamics is well defined 
whenever the initial distributions are stochastically dominated by some invariant measure 
v po . This simplifies also some technical estimates. One could remove the monotone as- 
sumption on g and choose other conditions assuring a well defined dynamics and some basic 
technical estimates involved in the proof, which would be similar to the ones appearing in 
[KL] [Chapter 5]. 

The entropy bound H{p N \v Pt ) ^ CoN d is rather restrictive. Indeed, given a locally 
Riemann integrable bounded profile po : R d — > [0, oo), let p N be the product measure on 
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p*j c M w ith slowly varying parameter associated to the profile po/m at scale N. Namely, 
p N is the product measure on N c((j) such that 

H N {r){x) = k) = Vp ( x /N)/m(v{x) = k ) ■ 

Due to the ergodicity of Q condition (|2.7p is fulfilled and, setting p' := sup x po(x), p N is 
stochastically dominated by v p ii m . On the other hand, the entropy H(p N \u Pt ) is given by 

w H(p M = w - [logy( - )-logy(^)j) + 

^ E {log^(p.))-lDgZ( V (^^)) (2.10) 

Hence, H(p N \p, Pt: ) ^ CoN d only if po approaches sufficiently fast the constant mp* at 
infinity. All these technical problems are due to the infinite space. In order to weaken the 
entropic assumption one should proceed as in |LMj . Since here we want to concentrate on 
the Moving Particle Lemma, which is the real new problem, we keep our assumptions. 

Finally, we have assumed uniqueness of the solution of differential equation (|2.9|) with 
initial condition po- Results on uniqueness can be found in |BCj . [KLJ [Chapter 5] and |Vaj . 
Proceeding as in |KL] [Section 5.7] and using the ideas developed below, one can prove 
that the limit points of the sequence {^Ite^T] are concentrated on paths t — > p(x,t)dx 
satisfying an energy estimate. 

3. Tightness of {w{*}te[o,T\ 

As already mentioned, in order to reduce the proof of Theorem 12. II to the Replacement 
Lemma one has to adapt the method of the corrected empirical measure developed in 
[GJ1| and after that invoke some homogenization properties proved in [F]. The discussion 
in |GJ1| refers to the zero range process on a finite toroidal grid and has to be modified 
in order to solve technical problems due to the presence of infinite particles. 

Given N £ N+, we define L^r as the generator of the random walk among random 
conductances to on the supercritical percolation cluster, after diffusive rescaling. More 
precisely, we define Cjv(w) = {x/N : x £ C(uj)} and set 

h N f(x/N) = N 2 £ u{x,x + e)[f{(x + e)/N)-f{x/N)} (3.1) 

eeB: 
x+e<=C(uj) 

for all x € C(oo) and / : Cn(lj) — > R. We denote by u~f the uniform measure on Cjv(w) 
given by = N~ d ^2 x ^c N (w) Below we will think of the operator L^v as acting on 
L 2 {v^). We write (•, -) v n and || • Hi^am for the scalar product and the norm in L 2 (v^), 
respectively. Note that L^r is a symmetric operator, such that (/, — L^r/)^ > for each 
nonzero function / G L 2 (v^). In particular, AI— L^v is invertible for each A > 0. Moreover, 
it holds 

(f,-h N g) u N = - E E + e )[f( x + e ) ~ f( x )] ■ [d( x + e ) ~ 9(x)] 

xeC(w) eeB: 

x+e£C(w) 
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for all functions f,g G L 2 (v*). Given A > 0, G G G£°(R d ) and JV G N+ we define G X N as 
the unique element of L 2 {v^) such that 

XG\ — LjvG^r = G A , (3-2) 

where G A is defined as the restriction to C N {oo) of the function AG - V • VVG G C^°(M d ). 
Let us collect some useful facts on the function G^: 

Lemma 3.1. Fix A > 0. Then, for each G G C™(R d ) it holds 

(G x N ,-h N G x N ) u M <c(G,A), (3.3) 

1 1 Gat 1 1 l 1 ^)) \\Gn\\l 2 (^) ^c(A,G), (3.4) 

ll^ivG^rllii^jv), ||LArG^|| i 2( J/ iv ) ^ c(A,G), (3.5) 

for a suitable positive constant c(A, G) depending on A and G, but not from N . Moreover, 
for Q-a.s. conductance fields uj it holds 

lim \\G X N - G\\» {vS) = 0, VG G C™(R d ) , (3.6) 
lim \\G X N - G\\ L , {uS) = 0, VG G C™(R d ) . (3.7) 

Proof. By taking the scalar product with G^ in (|3.2p one obtains that 



N 

M\Gn\\l 2 (vK) + (Gni -^nG^,) u n = (G^t,G A )^jv ^ ||G A / || L 2^jy)||G A || L 2 ( ^jy ) . 

Using that sup^y II C A II i 2 (^ t £^) < 00 > from the above expression one easily obtains the uni- 
form upper bounds on (Gjy, — LatG^)^ and HG^H^^jv). Since sup^ ^ i ||G a || L 2(^jv) < oo, 
by difference one obtains the uniform upper bound on HL/vG^H^ 



In order to estimate HG^H^i^jv) let us write Pt(x, y) for the probability that the random 
walk on Cn(uj) with generator L^r and starting point x is at site y at time t. Then, since 
the jump rates depend on the unoriented bonds, pf {x, y) = p^ (y, x). Since 

roo 

G X N (x)= / e- xt p?(x,y)G x (y), (3.8) 

y£C N (u)) 

for all x G Cn(u}), the above symmetry allows to conclude that 



i roo 

W G n\\li(»») < JfS E J e- xt p?(x,y)\G x (x)\ = 

£ |G A (x)|^i£|G A (n)|^<oo. (3.9) 

Again, since sup^f ^ x ||G A || L i( l/ jy) < oo, by difference one obtains the uniform upper bound 

on HLjvG^H^i^jv-j. 

The homogenization result (|3.6p follows from Theorem 2.4 (iii) in [F]. Finally, let us 
consider (j3.7|) . Given £ > 0, using Schwarz inequality, one can bound 



\Gn - GHii^iV) < 

\\G x N (u)I(\u\ > + \\G(u)I(\u\ > i)\\ L ^) + d d/2 \\G x N - G\\^ 
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Since G G C^°(M. d ) the second term in the r.h.s. is zero for I large enough. The last term 
in the r.h.s. goes to zero due to (|3.6p . In order to conclude we need to show that 



lim lim \\G x N {u)l{\u\ > t)\\ L i {vS) =0. (3.10) 



Since G x G C^(R d ) we can fix nonnegative functions F, / G C£°(R d ) such that -/ ^ G x ^ F. 
i\ f x 



We call Fh, the solutions in L 2 (v^) of the equations 



A-Fjy - = F . 

AT) ./jV 



respectively. From the integral representation (|3.8p we derive that F^, are nonnegative, 
and that —f^^ &n ^ ^jv on Cat (a;). In particular, in order to prove (|3.1Up it is enough 
to prove the same claim with F^, instead of G^. We give the proof for Fjy, the other 
case is completely similar. Let us define H G C°°(M. d ) as the unique solution in L 2 (dx) of 
the equation 

XH — V • V\7H = F . (3.11) 

Again, by a suitable integral representation, we get that H is nonnegative. Applying 
Schwarz inequality, we can estimate 

\\F x {u)l{\u\ > t)\W{v») = \\Fn\\lHu») - \\Fn(u)I(\u\ < £)\\i>(v») < 

pA II || L7/„,\TT/|„,| ^ 0\\\ I \\(TJ(~.\ IpA 



m»W) - \\H(u)I(\u\ < £)\\lMy3) + WiH{u) - F^)l{\u\ < e)\\ LHu N } ^ 
\\Fn\\v{u») - \\H(u)I(\u\ ^ t)\\LH^) + ci d l 2 \\F x - Hfg^ 



Since F^ and F are nonnegative functions, when repeating the steps in (j3.9j) with Ffo, F 
instead of G^, G x respectively, we get the the inequality is an equality and therefore 

II^/vIIl 1 ^) ^W^W^idx) = II-^IIl 1 ^) • 

This observation, the above bound (|3.12p and Theorem 2.4 (iii) in [F] imply that 

m\\F^(u)l(\u\>i)\\lA^) < \\H(u)I(\u\>e)\\iAW) 

At this point it is trivial to derive (|3,10p for F^. □ 

In the rest of this section, we will assume that w is a good conductance field, i.e. 
the infinite cluster C(ui) is well-defined and ui satisfies Lemma 13.11 We recall that these 
properties hold Q-a.s. 

The first step in proving the hydrodynamic limit consists in showing that the sequence of 
processes {Tr^jtgfo.T] i s tight in the Skohorod space D([0, T],M.). By adapting the proof of 
Proposition IV. 1.7 in |KLj to the vague convergence, one obtains that it is enough to show 
that the sequence of processes {7Tt V [G ! ]}tG[o,T] is tight in the Skohorod space D([0,T],R) 
for all G G C£°(R). A key relation between the zero range process and the random walk 
among random conductances is given by 

N 2 C (ir N (rj) [G]) = ± £ g( V (x)){lL N G)(x/N) . (3.13) 

The check of ()3.13p is trivial and based on integration by parts. At this point, due to the 
disorder given by the conductance field u, a second integration by parts as usually done for 
gradient systems (cf. |KL| [Chapter 5]) would be useless since the resulting object would 
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remain wild. A way to overcome this technical problem is given by the method of the 
corrected empirical measure: as explained below, the sequence of processes {^^[G]}^^^] 
behaves asymptotically as {^[G^JIjg^r] > thus the tightness of the former follows from 
the tightness of the latter. We need some care since the total number of particles can be 
infinite, hence it is not trivial that the process {^^ [G^]} t £[o y T] is we h defined. 
We start with a technical lemma which will be frequently used: 

Lemma 3.2. Let H be a nonnegative function on Cn{ui) belonging to Pi L 2 {v^) 

and let k ^ 0. Then 

sup W E Vt(x) k H(x/N) > a) < A- l c{k lPo )J\\H\\l 1{vS) + \\H\\ 2 L2(uS) 

(3.14) 

for all A > 0. 

Proof. We use a maximal inequality for reversible Markov processes due to Kipnis and 
Varadhan |KV] (cf. Theorem 11.1 in Appendix 1 of |KL]). Let us set 

F ^ = W E v(x) k H(x/N), (3.15) 

supposing first that H has bounded support. Note that F(rj) ^ F(rf) if rj(x) ^ rj'(x) for 
all x £ C(oj). Hence by the stochastic domination assumption, it is enough to prove (|3.14p 
with Pui,u P0 (always referred to the diffusively accelerated process) instead of P^w. We 
recall that u„ is reversible w.r.t. the the zero range process. Moreover 

v Po{F 2 ) = ^ Td Yl H{x/N)H{y/N)u p Mx) k n{y) k )^c{k,p )\\H\\ 2 L1{uSV (3.16) 
while 

TV 2 

Vpo (F, -N 2 CF) = J ^ £ H(x/N)H(y/N)v p Mx) k , -^(y) k ) < 

x,ydC(ui) 

N 2 

< k 'Po)^2d E E H(x/N)H(y/N). 

x£C{u)) yeC(u>): 
\x-y\=l 

Using the bound H{x/N)H{y/N) < H(x/N) 2 + H(y/N) 2 and the fact that d ^ 2, we 
conclude that 

u po (F,-N 2 CF) ^c(k,p )\\H\\ 2 L2{v „ y (3.17) 
By the result of Kipnis and Varadhan it holds 

P ^po( SU P Wd E Vt(x) k H(x/N) >A)^ ^tJv P0 {F 2 ) + Tu po (F, -N 2 CF) . 

(3.18) 

At this point the thesis follows from the above bounds (|3. 16[) and (|3.17p . In order to 
remove the assumption that H is local, it is enough to apply the result to the sequence of 
functions H n (x) := H(x)x{\x\ ^ n) and then apply the Monotone Convergence Theorem 
as n f oo. □ 
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Remark 1. We observe that the arguments used in the proof of Lemma 4-3 in [CLOj 
imply that, given a function H of bounded support and defining F as in (|3.15p , it holds 



sup (F( Vt ) - F(rj )) 2 ) ^ cTv po (F, -N 2 CF) 



s^T 



In particular, it holds 
E 



^ t ^ T 



sup F(th ? ^ cv po (F 2 ) + cTu po (F, -N 2 CF) . 



Using the bounds (|3.16p and (|3.17p . the domination assumption and the Monotone Con- 
vergence Theorem, under the same assumption of Lemma \3.S\ one obtains 



I sup 

■ ^ t ^ T 



1 



Wd Vt(x) k H(x/N) )^c(k,p ) 



(3.19) 



v o «: t t N a 



Using afterwards the Markov inequality, one concludes that 
Y, Vt(x) k H(x/N) > A) ^ c(k, Po )A- 

x<=C(ui) 

(3.20) 

for all A > 0. The use of (|3.14p or ()3.20p in £/ie resi of the discussion is completely 
equivalent. 

Due to Lemma [3.21 and Lemma [3.11 the process {^[G^]}^^] is well defined w.r.t 
K^^n. Let us explain why this process is a good approximation of the process {nf [G]}te[0,T\ : 

Lemma 3.3. Let G £ C c °°(]R d ). Then, given 5 > 0, it holds 



lim P, 



sup [G] \>S) =0. 

sg t ^ T 



(3.21) 



Proof. By Lemma 13.21 we can bound the above probability by 



N 



G\\ r1 /..«f\ + ||G^ 



At this point the thesis follows from Lemma I3TT 



N 



G\\ 2 



□ 



Due to the above Lemma, in order to prove the tightness of {7>"i V [G']} ig [o,T] it is enough 
to prove the tightness of [G^]}te[o,T} ■ Now we can go on with the standard method 
based on martingales and Aldous criterion for tightness (cf. [KLJ [Chapter 5]), but again 
we need to handle with care our objects due to the risk of explosion. We fix a good 
realization u of the conductance field. Due to Lemma I3.1|, Lemma I3.2I and the bound 
g(k) ^ g*k, we conclude that the process {M^}q < t < t where 

M t N (G):=K?(G x N )-Tr»(G x N )- J*^ g(r, s (x))h N G x N (x/N)ds, (3.22) 



is well defined P^iv-a.s. 
Lemma 3.4. Given 5 > 0, 

lim 



xec(ui) 



sup \M t N (G)\ ^ 5) = 



sC t < T 



(3.23) 
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Proof. Given n ^ 1, we define the cut-off function G^ n '■ CaK^) ~> K as G^ n (x) = 
G^{x)I{\x\ n). Then G x in is a local function and by the results of [X] (together with 
the stochastic domination assumption) we know that 

M^ n (G):=^(G^ n )-4{G^ n )- f± J2 g(Vs(x)% N G x N!n (x/N)ds 



N d 

xec(u}) 



is an L 2 -martingale of quadratic variation 

ft AT 2 



<M t N ' n (G)>=[ E 9 (ils(x)Mx,y)[G x N , n (y/N)-G x N<n (x/N)] 2 ds. 

xeC(ui) yeC(u>): 

\x-y\=l 

Note that, by the stochastic domination assumption and the bound g(k) ^ g*k, 

E^«(<M t N > n (G) >)< 
f* N 2 

9* N2dYl E E ^ P0 ^^)] w ( x ^)[ G iV,n(2//^)-^ in (x/iV)] 2 d S 

x£C(uj) yeC(u>): 

\x-y\=l 

= g* Po tN- d (Gi n , -1, n G x n „) u n . (3.24) 



By Doob's inequality and (|3.24p . we conclude that 

c 



sup \M t N ' n (G)\ > *) < >^(l<' n (G)| 2 ) < 



< t s£ T 

cg*T Po A A cVT>o ,„ 9t ^ 

Above we have used that lim nTOO (G^ n , — h^G^ n ) u N = (G^, — L,nGj^) u n < c(A) (see 
Lemma 13. ip . 

The above process {M t ,n (G)} tg [ 0j T] is a good approximation of {M® (G)} te \ 0j x] as 
n t oo. Indeed, it holds 



IfP^ sup \M?> n (G)-M»(G)\>6 =0, 5>0. (3.26) 

Indeed, since ||G^ n - G^IIl 1 ^)^ 2 ^) and ||L/vG^ n - JwvG^IIl 1 ^)^ 2 ^) conver g e to 
zero as N f oo and since ^ the above claim follows from Lemma 13.21 



At this point, fl52BD and (ET25]) imply (|3T23"1h □ 

Let us prove the tightness of {-k^ ' [G^Wt^^] using Aldous criterion (cf. Proposition 
1.2 and Proposition 1.6 in Section 4 of |KL| ): 

Lemma 3.5. Given G G C^°(W 1 ), the sequence of processes {^^[G^]}^^ ?>] ^ s tight in 
D([0,T],K). 

Proof. Fix 9 > and uppose that r is a stopping time w.r.t. the canonical filtration 
bounded by T. With some abuse of notation we write r + for the quantity min{r + #, T}. 
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Then, given e > 0, by Lemmata 13.11 and 13.21 

W E 9(r]s(x)% N G x N (x/N)d S 



xec(uj) 



>£ < 



(0<7* sup £ r7 a (x)|Lj V G^(x/JV)|tfa>e) <Cg*0/e. (3.27) 



In particular, 



lim lim sup sup P w „/ 

74-0 7v too Tj 6»g[0,7] 



^ E g(Vs(x))hNG x N (x/N)d S \> e) =0. (3.28 



An estimate similar to (|3.27p implies that 

^ 5 (r ?s (x))L^G^(x/iV)(i S 

x6C(w) 

Let us now come back to Lemma [3^41 Let r, as above. Then 



lim sup 

et°° AT > 1 



> e = 0. 



(3.29) 



limsupP^ u*r(\ M T+0 ~ M ? I > e )< lim sup „*r( sup |Mf| > e/2 1=0 (3.30) 

JVtoo ' ' jVfoo ' ^ s < T ' 



Collecting (|3.28[) . (|3.29p and (|3.30p . together with Lemma [331 we conclude that 



n? +e [G x N ]-ir?[G x )\>e)=0, 



lim lim sup sup 

7^0 Nfoo r,0G[O,7] 

lim sup r^ N {\nf[G x N }>e) =0. 

et°o AT > 1 / 

Aldous criterion for tightness allows to derive the thesis from (j3.31|) and (|3.32p . 



(3.31) 
(3.32) 

□ 



Let us come back to (|3.22p and investigate the integral term there. The following holds: 

Lemma 3.6. Let I(t) := /* N~ d ^ xeC{u}) g{rj s (x))(G^(x/N) - G{x/N))ds. Then, for all 
6 > 0. 

JfP„/( sup \I(t)\>6) =0. (3.33) 

Proof. Since g(k) ^ g*k and by Schwarz inequality we can bound 
m^J^Tg^Gb-Gf/* sup {N- d £ Va (x) 2 \G x N {x/N) - G{x/N)\}V 2 . 



< s < T 



x£C(lj) 



Using the stochastic domination assumption and applying Lemma 13.21 we obtain 



sup \I(t)\ >S) <W u , Vpo {J>8) < 



< t < T 



(l/*) 2 rV) 2 l|G$r - G\\ LH „N )y /\\G x N - G\\ 2 LHvS) + \\G X N - Gf~^ ) . (3.34) 



The thesis now follows by applying Lemma l3~7T 



□ 



ZERO RANGE PROCESSES ON THE SUPERCRITICAL PERCOLATION CLUSTER 



13 



We are finally arrived at the conclusion. Indeed, due to Lemma 13.31 and Lemma 13.51 we 
know that the sequence of processes {7iy }t e [o,T] is tight in the Skohorod space D([0, T],A4). 
Moreover, starting from the identity (|3.22p . applying Lemma l3~4"l using the identity (|3.2p 
which equivalent to 

h N G^ = XG X N -G x = \{G X N - G) + V • VVG , 

and finally invoking Lemma 13.61 we conclude that, fixed a good conductance field u, for 
any G G C™(R d ) and for any 5 > 

\G)-t$(G)- f ±2 ]T g( Vs (x))V-VVG(x/N)ds\>s)=0. 

(3.35) 

Using the stochastic domination assumption it is trivial to prove that any limit point of 
the sequence {tt^ }te[0,T] is concentrate on trajectories {^t}te[o,T] such that ir t is absolutely 
continuous w.r.t. to the Lebesgue measure. Moreover, in order to characterize the limit 
points as solution of the differential equation (|2.9p one would need non only (|3.35p . Indeed, 
it is necessary to prove that, given to good, for each function G G C^°([0, T] x R d ) it holds 



I™ F n sup 



HmV sup 7rf(G t )-<(G )- 

[ Jfd 9(Vs(x))V ■ VVG s (x/N)ds - J Trf (d s G s )ds\ > <f) = , (3.36) 

where G s (x) := G(s,x). One can easily recover (|3.36|) from the same estimates used to 
get (|3.35p and suitable approximations of G which are piecewise linear in t as in the final 
part of Section 3 in [GJ1J . In order to avoid heavy notation will continue the investigation 
of dSMD only. 

4. The Replacement Lemma 

As consequence of the discussion in the previous section, in order to prove the hydro- 
dynamical limit stated in Theorem 12.11 we only need to control the term 

I jfd E 9(Vs(x))V-VVG(x/N)ds. (4.1) 

To this aim we first introduce some notation. Given a family of parameters a%, oli, . . . , a n , 
we will write 

lim sup 

instead of 

lim sup lim sup • • • lim sup . 

Below, given x G and k G N, we write A x ^ for the box 

K,k := x +[-k,k] d DZ d , 
and we write r] k (x) for the density 

vk{x):= j2kTW £ v{y) ' 

yeA X)k nC(uj) 
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If x = we simply write instead of Ao,fc. 

Then, we claim that for Q— a.a. u), given G G C c (R d ), 5 > and a sequence \i 
of probability measures on N c ( w ) stochastically dominated by some u po and such that 
H(n N \v p .) ^ C N d , it holds 



limsupP^ n 



* 1 



yi/ Y g{Vs(x))G(x/N)ds- 



m N {x)/m)G{x/N)d S \ > s) = 0. (4.2) 



Let us first assume the above claim and explain how to conclude, supposing for simplicity 
of notation that eN G N. 

Given u G M. d and e > 0, define l UjE := (2e)~ d I{u G [— Then the integral 



^ l/z/W.e] > x ^ can De written as 



x/N,e\ 



i^ + lf N 
{2eNY V { h 



Let us define 



/ ^ E MK/NeMW -VVG{x/N)d8. 



(4.3) 



(4.4) 



Then, due to (|4.3p and since is Lipschitz with constant g*, we can estimate from above 
the difference between (|4.4p and the second integral term in (|4.2p with G substituted by 
V • VVG as 



(2eA^ + l) ' - (2eJV) £ 
(2eN) d 



f W E^)I V • Z?VG(*/J\0|cfe (4.5) 

" ™ 



Since the integral term in ()4.5p has finite expectation w.r.t 



and therefore also w.r.t. 



we conclude that the above difference goes to zero in probability w.r.t. 



At this point the conclusion of the proof of Theorem 12.11 can be obtained by the same 
arguments used in [KL| [pages 78,79]. 

Let us come back to our claim. Since 

Y g(Vs(x))G(x/N) = Y9(.Vb(x))I(x€C( U ))G(x/N), 

one can replace 



by a standard integration by parts argument and using that G G 
the first integral in (|4.2p by 

-t ^ % ^ ^ 

G{x/N) . 



IwT, [j2iTW ^ 9{Vs{y)) 

U x& d y ' y:yeA xJ nC(uj) 



Then the claim (14.21) follows from 



Lemma 4.1. (^Replacement LemmaJ For Q-a.a. oj, given 5 > 0, M G N and given 
a sequence of probability measures fx N on N c ^ stochastically dominated by some u po and 
such that H(^ N \u Pf ) ^ CoN d , it holds 



lim sup ' 

JVtoo,e4.0 



JO 



1 

~Nd 



Y V sN{r x 7] s ,T x uj)ds > S 



x-eA« 



0, 



(4.6) 
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where 



(2i + 1) 



y.y€AtnC(u) 



Let us define Tc 0j tv as the set of measurable functions / : N c ( w ) — > [0, oo) such that 
i) v pt (f) = 1, (ii) V(f) := Vp,(y/J,-Cy/J) < C N d - 2 and (iii) fdu Pt is stochastically 
dominated by u po (shortly, fdv Pt -< dv po ). Using the assumption H(p N \u Pt ) ^ CoN d and 
entropy production arguments as in [KL] [Chapter 5], in order to prove the Replacement 
Lemma it is enough to show that for Q-a.a. uj, given po, p*,Co, M > 0, it holds 

limsup sup / —7 V eN (T x r],T x u)f(ri)v pt (dr]) = 0. (4.7) 

Aftoo,e4-o/eT CoiJV J ^ „rr___ 



x£A M N 



Trivially, since v px stochastically dominates v p2 if pi > P2, it is enough to prove that, given 
P0) P*,Co, M > 0, for Q-a.a. uj (|4.7p is verified. We claim that the above result follows 
from the the One Block and the Two Blocks estimates: 

Lemma 4.2. (One block estimate,) Fix po, p*,Co, M > 0. Then, for Q-a.a. uj it holds 
limsup sup / V V e (T x r],T x uj)f(r])v pt (dri) =0. (4.8) 



x£A t 



Lemma 4.3. (Two blocks estimate) Fix po, p*, Co, M > 0. Then, for Q-a.a. uj it 
holds 

lim sup sup 

ATtoo,e|0/r°o/eT CoiJV 

[W E LeN + l¥ E \v l (y)-V £N (*)\)f(v> P Mv) = 0. (4.9) 



yeA 



We point out that the form of the Two Blocks Estimate is slightly weaker from the one 
in |KL] [Chapter 5] , nevertheless it is strong enough to imply, together with the One Block 
Estimate, equation (|4.7p . Indeed, let us define a(y) := I(y G C(uj)) and 



7d E \ Av yeA x , eN g(v(y))a(y) - A VyeA x , eN Av zeAy e g(ri(z))a(2 



N d 



Hv)--Jfd E |^eA*, e * (^eA^sM*))^) - m4>(r j i (y)/m) 



xGAmn 



J 3(^) := T^d E k%eA a: , eiV w^(r? l (y)/m) - m<j>{rf N (x) / m) 



N d 

xeA MN 

where denotes the standard average. Then 



(4.10) 
(4.11) 
(4.12) 



7d E V ^ij x r],r x Uj) 



N d 



xeAji 



^ \Av y£ ^ £N g{r ] (y))a{y)-m^ N {x)/m) ^ (h + I 2 + I 3 )(r,) . (4.13) 



xeAmn 
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Let us explain a simple bound that will be frequently used below, often without any 
mention. Consider a family of numbers b(x), x G Z d . Then, taking L,£ > we can write 

Av x£ a l H x ) ~ Av xe A L Av ue A x t b(u) = 



jl- &(»)(^6Ai)-r^rjt{«€A L :|!C-«| OQ <^}). 

XGA L+ £ 



In particular, it holds 



\Av xeAL b(x) - Av xeAL Av u£Ax A u )\ < rr-J E ( 414 ) 



,A L i 

xGA i+< .\A L _ f 

Due to the above bound we conclude that 

frfox jvd (2e ^ +1) d E E *>(«)«(«)• 

x€A MN u£A x ^ N+e \A XtSN _ e 

In particular, using that /<ifp t -< dz/ po , we conclude that J Ii{n)f{ji) v p*{drj) ^ c£/(eN) . 
The second term -^(t?) can be estimated for e ^ 1 as 



h(v)^j^ E A ^e^,ev v t(TzV>T z u) < jfi V e (T x r],T x uj) . 

x £A-ain 3;GA( M+1 ) JV 

Due to the One block estimate one gets that 

limsup sup / l2(rj)f(r])i>p t (dr]) = 0. 

7Vtoo,e4D/t°o /GT Co ,jv J 

The same result holds also for 13(77) due to the Lipschitz property of <p and the Two Blocks 
estimate. The above observations together with (|4.13p imply (|4,7D . 

5. Proof of the Two Blocks Estimate 

For simplicity of notation we set = (21 + l) d and we take M = 1 (the general case 
can be treated similarly). Moreover, given A C 7* d , we write M(A) for the number of 
particles in the region A, namely A/"(A) := YlxeAnCtu) r /( x )- 

Let us set 

■= -k E V.^h'M-fWl , (5-i) 



xGA 



.V 



B ^ := W E Av y^N A v zeAx , eN \v e (.y) ~ V e (z)\ , (5-2) 

x£A N 

C{1]) := w\a^\ E Av y^ x , eN E l^fo) - ^ w I • ( 5 - 3 ) 



|z-3/|oo > 2£ 

Since due to flUD) 

r? £jv (x) = A« 2GAa;iSjv ^(z)+5 , Ifl^clA^r 1 £ »/(«), 
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using that jdv Pf -< du po , in order to prove the Two Blocks estimate we only need to show 
that 

limsup sup / B(rf)f(r])up t (dr]) = 0. 

moo,eWtocf£Tc ,N J 

Since 

B{v) * c{v) + W\A^\ E AvyeA *-» E W(y) + v E (z)) , 

6 x£A N zeA X:EN : 

\z-y\aa < 21 

using again that fdv Pt -< dv po , in order to prove the Two Blocks Estimate we only need 
to show that 

lim sup sup 

JVtoo,e4.0,^too fer aQ ,N 

/f^wE E E K(V)-AA(A 2 ,)|/(r ? K(dr ?) =0. (5.4) 

|*-l/|oo > 2£ 

Let us now make an observation that will be frequently used below. Let X C Z d be a 
subset possibly depending on uj and on some parameters (for simplicity, we consider a 
real-value parameter L £ N). Suppose that for Q-a.a. w it holds 

nm sup —3 = U . 



Then, in order to prove (|5.4p we only need to show that 
lim sup sup 

/ wmW d ^ £ ^ K(a,,a^)-aa(a z ,a^)|/(^(^) = o. 

xeA N yEA x>eN zeA XjB n- 
\z-y\oo > 2< 

(5.5) 

We know that there exists ao > such that for each a G (0, ao] the random field Cj a 
defined in (|2.1|) is a supercritical Bernoulli bond percolation. Let us write C a (u>) for the 
associated infinite cluster. By ergodicity, 

limsup |(C \ C«) D Ajv|/|Ajv| =limQ(0 eC\C a ) = 0. (5.6) 

Hence, due to the above considerations, (|5.4p is proven if we show that, for each a 6 (0, ao], 
it holds (|5.4p with C replaced by C a . Moreover, since fdu Pt ~< v po , using Chebyshev 
inequality it is simple to prove that 

mJiwt^Atool N d (eN) 2d E E E 

x£A N yeA^sjsr z£A Xt£N : 
\z-y\ae ^ 2i 

I(Af(A y>i U K z/ ) > Aei)f( v )u Pt (drj) = . 
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At this point, we only need to prove the following: Fixed a G (0, ao] and A > 0, for Q-a.a. 
ui it holds 

\z-y\ao > 21 

\M(C a nA y/ )-M(C a nA z/ )\l(M(A y/ UA z/ )^A£i)f( v )u Pt (dr ] )=0, (5.7) 

where T^ o N is the family of measurable functions / : N c ^- 1 — > [0, oo) such that v Pt (/) = 1 
and V(f) k CoN d ~ 2 . 

We now use the results of [APJ about the chemical distance in the supercritical Bernoulli 
bond percolation ui a , for some fixed a G (0, ao]- We fix a positive integer K (this cor- 
responds to the parameter N in |APj . which is fixed large enough once for all). Given 
a E Z rf and s > 0, we set A as := A^jc+iWs- As in [APj . we call u a the microscopic ran- 
dom field. The macroscopic one a = {cr(a) : a G Z d } G {0, is defined in |APj stating 
that cr(a) = 1 if and only if the microscopic field w Q satisfies certain geometric properties 
inside the box A a ^ K / A . These properties are described on page 1038 in |APj . but their 
content is not relevant here, hence we do not recall them. What is relevant for us is that 
there exists a function p : N — > [0, 1) with lim^ 0O p(ir) = 1, such that a stochastically 
dominates a Bernoulli site percolation of parameter p(K) (see Proposition 2.1 in [APj ) . 
Below we denote by Wp(K) the law of this last random field, taking K large enough such 
that p(K) is supercritical. As in [AP] we call a point a G 7j d white or black if <x(a) = 1 or 
respectively, and we write in boldface the sites referred to the macroscopic field. Recall 
that a subset of Z d is *-connected if it is connected with respect to the adjacency relation 

x ~ \x - y\oo = 1 . 

C* is defined as the set of all *-connected macroscopic black clusters. Given a G C* 
denotes the element of C* containing a (with the convention that C* = if a is white), 
while C* is defined as C* = C* U d out C* a . We recall that, given a finite subset A C Z d , 
its outer boundary is defined as 

d out A := {x G A c : 3y G A, {x, y} G E d } . 

We use the convention that d out C a = {a} for a white site a G Z, d . Hence, for a white it 
holds C* = {a}. 

Let us recall the first part of Proposition 3.1 in [AP]. To this aim, given x,y G 
we write a(x) and &(y) for the unique sites in Z d such that x G A ai x and y G A a> ^. 
We set n := |a(x) — a(y)|i and choose a macroscopic path A XjV = (ao,ai, . . . ,a n ) with 
ao = a(x) and a n = a(y) (in particular, we require that |aj — aj + i|oo = 1). We build the 
path A Xj y in the following way: we start in a(x), then we move by unitary steps along 
the line a(x) + Zei until reaching the point a' having the same first coordinate as a(y), 
then we move by unitary steps along the line a' + Ze2 until reaching the point having the 
same first two coordinates as a(y) and so on. Then, Proposition 3.1 in [AP] implies (for 
K large enough, as we assume) that given any points x,y G C Q there exists a path ^ x ^ y 
joining x to y inside C a such that j x>y is contained in 

W X)V := UagA^ (u we c.,A Wi5 ^ /4 ) . (5.8) 
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These are the main results of [APj that we will use below. Note that, since the sets C* 
can be arbitrarily large, the information that ^ X)V C W x y is not strong enough to allow to 
repeat the usual arguments in order to prove the Moving Particle Lemma, and therefore 
the Two Blocks Estimate. Hence, one needs some new ideas, that now we present. 

First, we isolate a set of bad points as follows. We fix a parameter L > and we define 
the subsets B(L),B(L) C Z d as 

B(L) := {a G Z d : |C*| > L} , (5.9) 
B{L) :=U aeB (L)A a ,io^. (5.10) 

Lemma 5.1. Given a in (0, cto], for Q-a.a. uj it holds 

v \B(L) HA N \ 

hmsup — — : = 0. (5-11) 

Aftoo,Lt°o |Ajv| 

Proof. Since \B(L) n Ajy\ ^ c(-?T)|B(L) n Ajv|, we only need to prove the thesis with B(L) 
replaced by B(L). We introduce the nondecreasing function pi : N — > [0, oo) defined as 
Pl{p>) := > L)n. Then we can bound 



|B(L) n Ajv| Y, Pl( C *) 



C*nAjv^0 

Since a stochastically dominates the Bernoulli site percolation with law 'Pp(K) an d due to 
Lemma 2.3 in |DPj . we conclude that 



B(L)nA N \ >a\A N \) < q( Y p L (C*) > a\A N \) ^ 



CnAjv^0 

- rtK<[ J2 PL(C*)>a\A N \) <P( PL(Cl)^a\A N \), (5.12) 

C*GC*: aeAjv 

where the random variables C* (called pre-clusters) are i.i.d. and have the same law of 

Cq under Pp(^). Their construction is due to Fontes and Newman |FNlj . |FN2j . Due to 
formula (4.47) of |AP| . £f>_ (K) (|C5|) is finite for K large, in particular 

YimE(p L {C* )) = Q. (5.13) 

By applying Cramer's theorem, we deduce that 

P( Y, PLi&L) > 2E( PL (& ))\A N \) < e' c ^ Nd , 



N 



for some positive constant c(L) and for all N ^ 1. Hence, due to (|5.12p and Borel-Cantelli 
lemma, we can conclude that for Q-a.a. uj it holds 

\B(L) n A N \/\A N \ < 2E( PL (C* )) , ViV > N (L,u) . 

At this point, the thesis follows from (j5. 13|) . □ 
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At this point, due to the arguments leading to (|5.5p . we only need to prove the following: 
given a G (0, oto] and A > 0, for Q-a.a. to it holds 



h ™l u p ™? / N d(Xf\2dpd E E E 



7Vtoo, e wt°o,Ltoc /eTf. „ 7 N d (eN) 2d £ d 

\z-V\oa > 2t 

W(T yAa ) - M(T zAa )\l(M(T yAa U T zAa ) ^ A£i)f( V )u p , (dr,) = (5.14) 

where 

T uAa = (A u/ nC a )\B(L), u£Z d . (5.15) 

Above we have used also that 

I(Af(A yi£ U A z/ ) ^ Alt) < I(Af(T yAa U r^, a ) < Alt) ■ 

Note that in the integral of (|5,14p . the function / multiplies an J^-measurable function, 
where Fn is the c-algebra generated by the random variables {rj(x) : x G Gn} and Gn 
is the set of good points define as 

G N := (A N+1 DC a )\B(L). (5.16) 

Since P(-) is a convex functional (see Corollary 10.3 in Appendix 1 of |KLJ), it must be 

V{u p ,(f\^N)) <£>(/) ^C N d - 2 . 

Hence, by taking the conditional expectation w.r.t. J~n in (j5. 14|) . we conclude that we 
only need to prove (|5.14|) by substituting Y^ o N with Tjn. n defined as the family of 
measurable functions / : N c(w) -)• [0,oo) such that u(f) = 1 and V(f) < C N d - 2 . 

Recall the definition of the function ip(-) given before (|2.4p . By the change of variable 
r\ — > n — 5 X one easily proves the identity 

'9(v*)(Vf(rf l > y ) - \/7(r?)) 2 ] = v{p*)v P « [{x/f(v x ' + ) - x/f(v y ' + )) 2 } , (5.17) 

where in general n z,+ denotes the configuration obtained from rj by adding a particle at 
site z, i.e. n z ' + = r/ + 5 Z . Let us write V x , y for the operator 

V.Xr?) := h{^ + ) - h(r?>>+) . 

We can finally state our weak version of the Moving Particle Lemma: 

Lemma 5.2. For Q-a.a. uj the following holds. Fixed a G (0, cto] an d L > 0, there exists 
a positive constant k = k(L, a) such that 



W* 



E E E E E 



N d {eN) 2d £ 2d * 

\Z-y\ao > 21 

^N 2 - d V(f)/K^C /K, (5.18) 

for any function f G YjL n and for any N, £, Cq. 

Proof. Recall the definition of the path j Xjy given for x,y G C a in the discussion before 
(|5.8p . Given a bond b non intersecting Gn, since / is J-jv-measurable it holds Vb\ff = 0. 
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Using this simple observation, by a standard telescoping argument together with Schwarz 
inequality, we obtain that 

n— 1 n— 1 

Vp. ((V^a//) 2 ) < { K{ui,u i+ i} n G N + 0)} • { v P * ((V« ilUi+1 \/7) 2 ) } , (5.19) 

i=0 i=0 

where the path ^y U:V is written as (it = uo, u\, . . . , u n = v). Recall that if v ptr ((V Uz ,u i+1 Vf) 2 ) ^ 
then the set {ui,Ui + i} must intersect the set of good points G/v defined in (|5. 16j) . 

If b is a bond of j u ,v, then b must be contained in the set W UfV defined in (j5.8j) . In 
particular, there exists a G A UjV and w G C* such that b is contained in A Wi5 ^ / / 4 . Denoting 
doo(-,-) the distance between subsets of 7L d induced by the uniform norm | • | ^ , we can 
write 

d^b, {2K + l)w) 5K/4. (5.20) 

We claim that, if b intersects Gat, then |C*| ^ L. If a is white then |C*| = 1 and the 
claim is trivally true. Let us suppose that a is black and that |C*| > L. By definition of 
C*, there exists some point a' G C* such that |w — a'|oo ^ 1, i.e. 

doo {{2K + l)w, (2K + l)a') 2K + 1 . (5.21) 

Due to (|5.20p and (|5.2ip we conclude that b C A^io^- = A-(2K+i)a>,ioK- On the other 
hand, C*, = C* and by definition of the set of bad points we get that A a / t iQK C B(L). The 
above observations imply that b C B(L) in contradiction with the fact that b intersects 
the set of good points Gjy. This concludes the proof of our claim: |C*| ^ L. 

We define |7«,u|* := Y^l=o ^-{{ u ii u i+i} H Gn ^ 0)- We claim that for almost all conduc- 
tance field uj there exists a constant c(K, L) depending only on K and L such that, for all 
(x,y, z,u,v) as in (j5. 18|) . |7u,u|* ^ ceN . Indeed, by the above claim, we get that 

\luA* ^ |Asj?/4| l^aWlC*! ^ L) ^ c{d)L\K bK/A \ \A UiV \ ^ C(K,L)eN. 

Given a bond &SE^ let us estimate the cardinality of the set X(6), given by the strings 
(x,y, z,u,v) with x,y,z,u,v as in the l.h.s. of (j5. 18fl . such that b is a bond of the path 
7 Mjl) and b intersects G^. Up to now we know that there exist w,a such that b intersects 
A W)5 #y4, a G A u ^ v , w G C* and |C*| ^ L. In particular, it must be 

doo(6, {2K + \)A U>V ) ^ 

doo(6, (2K + l)w) + (2K + 1)^ (w, C a ) + (2K + l)diam(C;) ^ 

hK/A + {2K + !){!+ L). (5.22) 

Hence, if (x,y,z,u,v) G X(6), then the distance between b and {2K + l)A UjV is bounded 
by some constant depending only on K and L. Note that the macroscopic path A u>v 
has length bounded by ceN/K. Let us consider now the set Y(b) of macroscopic path 
(ao, ai, . . . , a&) such that 

k s: ceN/K , (5.23) 

d 00 (b,(2K + l){a , ai ,...,a k }) ^ 5K/A + (2K + 1)(1 + L) . (5.24) 

By the same computations used in the proof of the standard Moving Particle Lemma, one 
easily obtains that |Y(6)| ^ c(K, L)(eN) d+1 . Fixed a macroscopic path in Y(6) there are 
at most c(K) ways to choose microscopic points u, v such that the macroscopic path equals 
A UjV . Fixed also u and v, we have at most i"l d ways to choose (y, z) as in (|5.18p . Fixed (y, z) 
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we have at most c(eN) d ways to choose x as in ([5TT8]) . Hence |X(6)| < c(K, L)(sN) 2d+1 £ 2d . 
Finally, recall that the length of 7 Uj „ is bounded by ceN . Since moreover all paths 'ju.v 
are in C a , from the above observations and from (|5. 17|) we derive that 

i.h., ofvm < E mMMVD 2 ) < 

r(K T 1 \(rN) 2d + 2 r- 2 P 2d 

N d (sN) 2d i 2d a * V{f) * C{K ' L ' a ^ dV ^ < ^ L ' ' ( 5 ' 25 ) 

□ 

Since a is fixed, we will stress only the dependence of the constant n(L, a) in Lemma 
on L writing simply k(L). We introduce the function 

Z< ti) := Nd( £ N\2d EE E 

xeA N yeAx te n zeA x>eN : 
\z-y\oo > 2£ 

Af(T yAa ) Af(T zAa ) 



I(M(T yAa U T zAa ) < A^) . 

Recall that we need to prove (|5.14[> where T^ Q N is substituted by the family Tfjl, N 
defined before (|5.17p . Due to the above Lemma, we can substitute TjL n with the family 
C /k(l) n °f measurable functions / : N c ( w ) — > [0, oo) such that v{f) = 1 and such 
that the l.h.s. of (|5.18p is bounded by Cq/k{L). Namely, we need to prove that, given 
a S (0, ckq] and A > 0, for Q-a.a. u it holds 



limsup sup u pt (Z(j])f(r])) = . (5.26) 

JVt« 

Since by Lemma 15.21 



NI-oo.elOltoo.Ltoo f(z -r^ 

/fcl C / K (L),iV 



sup ^ P AVf,- N dr £ N)2d £ 2d E E E E E v ^\/7)^^, 

/ eT c /K(i),iv * ^eAjv 2/6A x , e jv zeA^^jv: uer Si « iQ , uer^c 

we only need to prove that, given a £ (0, cto] and A, 7 > 0, for Q-a.a. uj it holds 

limsup sup spec L 2( u Jz(r]) + V V V u ,v } < , (5.27) 

where sup spec L 2r v \{-) denotes the supremum of the spectrum in L 2 {v Pt ) of the given 
operator. Now we use the subadditivity property 
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to bound the l.h.s. of ([OZJ) by 

JV(r«j,a) AT(r., a ) 



^tooSJS^toolvW E E E -p^(^){ 

\z-y\ x > 21 

7 «(L) 



/>d /d 



W, Aa ur^)<^ + |U £ £ V „ )t ,}. (5.28) 

We observe that the operator inside the {-}-brackets depends only on r/ restricted to T y>z := 
r y /, Q ur 2j £ iQ . By calling Vk,y,z the canonical measure on <Sfc,y,z '■ = 

obtained by conditioning the marginal of v Pit on N Ty ' z to the event {M(T ytZ ) = k}, we can 
bound (|5lZ8"j) by 



E E E 



|2-S/|oo 3* 2£» 

supspec L2(lw) ( ^ ^ +72^" E E V ^j- ( 5 - 29 ) 

Given integers k,ni,ri2 G N, define for j = 1,2 the set := {1, 2, . . . , n^}, with the 
convention that Tj = if rij = 0. Then define the space 

S k , ni ,n 2 ■= {(Ci,C 2 ) G N Fl x N r2 : £ Ci(oi) + E ^(« 2 ) = k ^ 

aieVi a 2 er 2 

and set J\f(d) = X^er^ Ci( a i)- Finally, call 

Vk,n\,ri2 the probability measure on Sk,m,n 2 
obtained by first taking the product measure on N Fl x N r2 with the same marginals as v pit , 
and afterwards by conditioning this product measure to the event that the total number 
of particles is k. Finally, define 

F{k, ni ,m,e,£,L) :=supspec i2( ^ iniiii2) {c d K(Ci)-AA(C 2 )| + ^§ Yl E V «>4 • 

* «eri ^er 2 

(5.30) 

Note that the operator r yn(L)l~ 2d ^ ugri S t) £r 2 ^u,v is the Markov generator of a process 
on Sk : m,n 2 such that the measure i>k,n±,n 2 is reversible and ergodic. In particular, is a 
simple eigenvalue for this process. Fixed £, we will vary the triple (k, m, n 2 ) in a finite set, 
more precisely we will take nx,n 2 ^ ^* and ^ k ^ A£f.. Then, applying Perturbation 
Theory (see Corollary 1.2 in Appendix 3.1 of |KLj ). we conclude that 



lim sup sup 

e^O k,ni,n 2 



F(k,rii,n 2 ,s,£,L) — G(k,£,ni,n 2 ) =0 (5.31 



where 

G(k,£, ni ,n 2 ) = v K n u n 2 {l* d \M{ti) - N(Q 2 )\) ■ (5.32) 

The above result implies that in order to prove that (|5.29p is nonnegative we only need to 
show that 

limsu P xrdftxrW E E E sup G(k,£,n y ,n z ) < 0, (5.33) 
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where 



fly — l^y^cul , Tl z — |r 2i £ )Q ,| . 



Lemma 5.3. Given 5 > 0, for Q-a.a. u it holds 

l{\n y /ii -m a \>5 or \n g /e* - m a \ > 8) = , (5.34) 

where m a = Q(0 G C a ). 

Proof. Recall definition ()5.15p and set N y = \C a fl Ay ^|, = |C a n A 2j ^|. Then 

]^vf E E E i(k-^i>^)^ 



x&A N j/GA^^jv zeA^.ejv 

c 



N d (eN) c 



xeAjv yGA 2 



E E l^)nVI%WnA 



N+eN\ 



5N d (eN) d £ d V 1 v ' W|x ^iV^ 1 

Due to the above bound, the same bound for n z and N z , Lemma 15. II and finally a 5/2- 
argument, we conclude that it is enough to prove (|5.34p substituting n y and n z with N y 
and N z respectively. Moreover, due to ergodicity, for Q-a.a. ui it holds 



limsup Ndflmd E E K\N y /£i-m a \>6) 



xgAjv yeA x>eN 



limsup — , V I(|JV„/^ - m a | ><5) = c limsup Q(|Ay^ - m a \ > 5) =0. 



3/GA 2 jv 



£foo 



A similar bound can be obtained for z instead of y, thus implying (|5.34p . 
Since G(k,£, n\, n 2 ) ^ A and since 

AA(Ci) k 711 



(5.35) 
□ 



G(k,t,m,n 2 ) < ffc,ni,n 2 ( 



m + n 2 

' AA(C 2 



+ 



k n 2 



m + n 2 



+ ^^^, (5.36) 
ni + re 2 



by the above Lemma in order to prove (|5.33p we only need to prove for j = 1,2 that 



limsup sup Vk,m,n 2 [ 

et°o,8l0 (k,m,n 2 )eX V 



k rii 



n y + n 2 i d 



0. 



(5.37) 



where 
X = {(fc,ni,n2) G N 3 : 

|j < i4 , |e[m a -5,(m Q + ^Al], ^ e [m a - 5, (m Q + 6) A 1]} . (5.38) 

At this point (|5.33p derives from the local central limit theorem as in in [KL] [page 89, 
Step 6]. 
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6. Proof of the One Block Estimate 



We use here several arguments developed in the previous section. In order to avoid 
repetitions, we will only sketch the proof. As before, for simplicity of notation we take 
M = 1. 

Let us define m a := Q(0 G C a ). Note that m = lim^o 77 ^- Setting a = n e (y), 
b = M(C a n A y ^)/£ d , since <fi is Lipschitz we can bound 



m(j)( — ) — m a (f)( 

m m a 



^ \m - m a \(j){—) + m a \4>(—) - 

m mm 



— )l + m a \0(— ) 



m — m, 



Ct I * * ^^(X I 7 i J 

— g a + g \a — o\ + g bm a 

m m 



1 1 

m m a 



\m-m a \ „ t 



m 



g*rf{y) + g*^M(A y/ \ C a ) + g*M(C a n A y/ )£7 d m a 



m 



m m c 
1 1 



m 



m m r 



Note that Q{rf) is increasing in r\. Hence, using that f(r])vp t (dr]) -< f(n)^ po (dr]), we easily 
obtain that 

limsup sup I — ^ V g(T x r/,T x uj)f(r])u pt (dr]) =0. (6.1) 
Ntoo,etooMOfeT Co ,N J iV ST 



Due to the above result and reasoning as in the derivation of (|5.14p where Y£, Q N can 

be replaced by T|L o n (see the discussion after (|5.14p ). we only need to prove that given 
C > 0, A > and a G (0, a ], for Q-a.a. oj it holds 

lim sup sup 

Av x£ a n 

C d E 9Wv))-rrhx<f>(M(Txte)/mc,t§ \l(M(r x4 , a ) < Aif) f ( V )u^( V ) = 0. 

(6.2) 

At this point by the same arguments of Lemma 15.21 one can prove that 

Lemma 6.1. For Q-a.a. uj the following holds. Fixed a G (0, cto] an d L > 0, there exists 
a positive constant k = k(«, L) such that 

p(p*)N 2 l: 2d - 2 Av xeAN E *V(( V ^'V7) 2 ) < N- d+2 V(f)/K^C /K, (6.3) 

for any function f G YjL N and for any N, £, Cq. 

Due to the above lemma, as in the derivation of (|5,27p we only need to prove that given 
a G (0,ao], A, j > 0, for Q-a.a. uj it holds: 

lim sup sup spec L 2 u, \ 

N^oo,t[oo,L"[oo 

{av x&An |c d Yl 9(v(y)) - m a <f>(jV(r xAa )/mjf) \i{M(r xAa ) < Alt) 



+ ^ E E V„»]}<0. (6.4) 
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Using subadditivity as in the derivation of (|5.28p we can bound the above l.h.s. by 



limsup Av xe A N sup spec L 2( Vp ^ 



{ |C d E 9{V(V)) - m a ^(N(T xAa )/m a li) \l{M{T x ^ a ) ^ Aif) 



s/erv, 



^-W- E E V„]}- (6.5) 



Again by conditioning on the number of particles in T* x g a and afterwards applying per- 
turbation theory (see (|5.29p . (|5.3ip and (|5.33p ). one only needs to show that 

(6.6) 



limsup Av x( z\ N sup G(k,£,n 1 

Aftoo,^too,Ltoo fce{0,l,...,A^} 



0. 



where n x := |r X)Qj £|, v^ n is the measure on {( G N" : Ya=i CW = ^} obtained by taking 
the product measure with the same marginals as u pt and then conditioning on the event 
that the total number of particles equals k, and where 

n 

G(k, I, n) := u Kn [|C d ~ m a 0(fc / 'mj*) 



i=l 



One can prove that for Q— a.a. u it holds 

limsup Av xeAN K\n x /ii - m a \ > 5) = 0, 

for each positive constant 5 > 0. As in [KL] [Chapter 5], one has that 



lim sup i/fc n 



n 



4 EfK«)-#/")l 



i=l 



(6.7) 



(6.8) 



where 



J:= |(fc,n) eN 2 — a,^£ (m a - <5, (m a + 5) A 1)} 
At this point (|6.6|) follows from (|6.7|) and 



Appendix A. Zero range process on Z d with random conductances 

Recall that the enviroment uj = (u(b) : b G Ed) is given by a family of i.i.d. random 
variables parameterized by the set Ed of non-orientied bonds in Z d , d ^ 2. We denote by 
Q the law of uj, we assume that Q(w(6) 6 [0, co]) = 1 and that Q(u(b) > 0) is supercritical. 
We fix a function g : N — > [0, oo) as in Subsection 12.21 Given a realization of uj, we consider 
the zero range process rj(t) on Z rf whose Markov generator N 2 C acts on local functions as 

N 2 Cf( V ) = iV 2 E E sM^M*. x + e) {f( V x ' x+e ) - f(r,)) , (A.l) 

e€B x£Z d 

where 5 = {±ei, ±e2, . . . , ±ed}, e\, being the canonical basis of Given an 

admissible initial distribution p, N on {0, l} z<i (i.e. such that the corresponding zero range 
process is well defined), we denote by P w njv the law of : i ^ 0). Trivially, the zero 
range process behaves independently on the different clusters of the conductance field. 
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If Q(w(6) > 0) = 1, then Q-a.s. the infinite cluster C(uj) coincides with and the 
hydrodynamic behavior of the zero range process on Z d is described by Theorem 12.11 If 
Q(uj(b) > 0) < 1 the bulk behavior of the zero range process on Z d is different due to the 
presence of finite clusters acting as traps, as we now explain. 

First, we observe that the finite clusters cannot be too big. Indeed, as byproduct 
of Borel-Cantelli Lemma and Theorems (8.18) and (8.21) in [G], there exists a positive 
constant 7 > such that for Q-a.a. uj the following property (PI) holds: 
(PI) for each N ^ 1 and each finite cluster C intersecting the box [— N, N] d , the diameter 
of C is bounded by 7ln(l + N). 

Lemma A.l. Suppose that uj has a unique infinite cluster C{uj) and that uj satisfies the 
above property (PI)- Let G G C c (M. d ) and 77 G N zd be such that n(x) = for all x G C(oj). 
Call Ac the support of G and call Aq the set of points z £ M. d having distance from Aq 
at most 1. Then there exist positive constants iVo(G, 7), C(G, 7) depending only on G 
and 7 such that the zero range process on Z d with initial configuration r\ satisfies a.s. the 
following properties: r]t(x) = for all x G C(uj) and, for N ^ Nq, 

N- d £ G(x/N) V (z) - N- d £ G(x/N) Vt (x) < C(G, 7 )^^ £ r,(x) . 

x&L d x&L d " x£Z d :x/NeA G 

(A.2) 

We point out that the zero range process is well defined when starting in r), indeed the 
dynamics reduces to a family of independent zero range processes on the finite clusters, 
while the infinite cluster C(oj) remains empty. 

Proof. The fact that rjt(x) = with x G" C(uj) is trivial. Let us prove (fO]) . Without loss 
of generality we can suppose that G has support in [— 1,1] (the general case is treated 
similarly). Let us write C± , C^, . . . , C^ N for the family of finite clusters intersecting the 
box [— N, N] d . For each cluster C^ we fix a point xf G . Since by the property (PI) 
each Cf has diameter at most 7ln(l + N), we have 

\G(x/N) - G{xf /N)\ < C(G)jln(l + N)/N , Vi : 1 < i< k N , Vx G Cf . 

The above estimate implies that 



N- d Y,G{x/N) m (x)-N- d Y,G{x?/N) £ Vt (x) 
x&z d i=l xeC^ 



k N 

t=l iEgCf 

Using now that the number of particles in each cluster is time-independent and that for 
N large enough Cf C [— 2JV, 2N\ for alH = 1, . . . , kj^, we get the thesis. □ 

As a consequence of Lemma lA.li if uj has a unique infinite cluster and if uj satisfies 
property (PI), then for any admissible initial configuration tjq (i.e. such that the zero 
range process on Z d is well defined when starting in 170) and any G G G c (R d ), it holds 

N- d G(x/N)n t (x) = N- d £ G(x/N)r] t (x) + N~ d £ G(x/N) m (x) + o(l) , 

(A.3) 
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where |o(l)| ^ C(G, 7) ln ^ d +f } Y,x&iA:x/N&K G 7 ?( x )- At this point, denoting by fj, N the 
initial distribution of the zero range process r\t on one can derive the hydrodynamic 
limit of r\t if the marginals of p, N on C(oj) and Z d \ C(to), respectively are associated to 
suitable macroscopic profiles. In what follows, we discuss a special case where this last 
property is satisfied. 

We fix a smooth, bounded nonnegative function po : M. d — > [0, 00) and for each N we 
define £t N as the product probability measure on N zd such that for all x G Z d it holds 

fl N (r](x) = k) = v Pq ( x /n){v(x) = k) , (A.4) 



where v p is defined as in Subsection 12.21 with the difference that now it is referred to all 
Z d and not only to C{to). 

We call p N the marginal of p N on C(oj): p N is a product probability measure on 
satisfying (|A.4j) for all x G C(w) (note that p^ depends on w). Similarly we call ^ Pi c(uj) 
the marginal of v p on C(w). By the discussion at the end of Section [21 if the smooth profile 
po converges sufficiently fast at infinity to a positive constant p*, it holds 

limsupN- d H(p N \v pt C{uj) ) < 00 Q-a.s. (A. 5) 

Vtoo 

Theorem A. 2. Suppose that the bounded smooth profile po : M. d — > [0, 00) satisfies (|A.5[) . 
Then for all t > 0, G G C c (M d ) and 5 > 0, for Q-a.a. ui it holds 

N- d V G(x/N)th(x) - [ G(x)p(x,t)dx > S) = (A.6) 
where, setting m = Q(0 G 

p(x, i) = mp(x, t) + (1 — m)p {x) (A. 7) 

and p : x [0, 00) -> 1 is the unique weak solution of the heat equation 

8 t p = V- (VV<f>(p)) (A.8) 
with boundary condition po = p ai i = 0. 

Proof. Since p, N is a product measure and the dynamics on different clusters is indepen- 
dent, the process restricted to C(to) has law P^w. As discussed at the end of Section [21 
to this last process we can apply Theorem 12. II Since, for Q-a.a. to the initial distributions 
p N are associated to the macroscopic profile mpo, we conclude that for Q-a.a. to it holds 



Nfoo ^ 



V G(x/N)rj t (x)-m f G(x)p(x,t)dx >s)=0. (A.9) 



Let us now consider the evolution outside the infinite cluster. Let us write 

N- d G(x/N) m {x) = N- d Y J G{x/N) m (x)-N- d £ G{x / N) ^{x) . 

x£C(uj) x£Z d xGC(uj) 

We know that, when n is sampled with distribution p, , the addenda in the r.h.s. converge 
in probability to J G(x)po(x) and m f G(x)po(dx), for Q-a.a. to. As a consequence the 
l.h.s. converges in probability to (1 — m) J G(x)po(x) for Q-a.a. to. In addition, 



sup { / p N {dr,)N- d V n(x)\ <oo. 

N^l ^ J f-f w > 

x-.x/NeAo 
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The above observations and Lemma lA. II (cf. ()A.3j) ) imply that 



lim P. 

iVfoo 




The thesis then follows from (|A~9|) and (fATTOj) . □ 
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